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We present a microscopic theory of the long-wavelength magnetic properties of the ferromagnetic 
diluted magnetic semiconductor (Ga,Mn)As. Details of the host semiconductor band structure, de- 
scribed by a six-band Kohn-Luttinger Hamiltonian, are taken into account. We relate our quantum- 
mechanical calculation to the classical micromagnetic energy functional and determine anisotropy 
energies and exchange constants. We find that the exchange constant is substantially enhanced 
compared to the case of a parabolic heavy-hole-band model. 

PACS numbers: 75.50.Pp, 75.30.Ds, 75.30.Gw 



I. INTRODUCTION 

The recent discovery of carrier-mediated ferromag- 
nctism in diluted magnetic semiconductors (DMS) has 
generated intense interest, in part because it suggests 
the prospect of developing devices which combine in- 
formatioB, psocessing and storage functionalities in one 
material .ErEil Ferromagnetism has been observed in Mn 
doped GaAs up to critical temperatures Tc of llOK (see 
Ref. |ll|). Doping a III-V compound semiconductor with 
Mn introduces both local magnetic moments, with con- 
centration NyY-a and spin S = 5/2, and itinerant valence- 
band carriers with density p and spin s = 1/2. An an- 
tiferromagnetic interaction between both kinds of spin 
mediates an effective ferromagnetic interaction between 
Mn^^ spins. 

A phenomenological long-wavelength description of 
ferromagnets usually requires only a small set of char- 
acteristic parameters. For example, a ferromagnet which 
possesses an uniaxial anisotropy can be modeled by a 
classical micromagnetic energy functional £^[n(r)]. 



£;[n(r)] = 



Ksiii^e + A(yhy 



(1) 



where n(r) is the unit vector that specifies the (space- 
dependent) local Mn spin orientation, and 9 is its angle 
with respect to the easy axis. The dependence on the ori- 
entation of the magnetic moment is parametrized by the 
anisotropy constant K. The exchange constant (or spin 
stiffness) A governs the energy cost to twist the orienta- 
tions of adjacent spins relative to each other. Together 
with the saturation moment fJ-oMg and the magnetostatic 
energy, dropped for convenience in Eq. (|^), the param- 
eters K and A determine a whole variety of magnetic 
properties such asc3 domain- wall width 6b , domain- wall 
energy per area 7, exchange length Zex, hardness param- 
eter K, single-domain radius Rsd, and anisotropy field 
HoHq. In addition they determine the energy cost fl^ of 



collective long-wavelength spin excitations, spin waves. 
For uniaxial anisotropy the quantized energy is 



2K 



2A 



(2) 



at wavevector k. Here, g is the g-factor and fiB the Bohr 
magneton. There is an energy gap which is determined by 
the anisotropy K. The spin exchange constant A defines 
the curvature of the spin-wave dispersion. 

The purpose of this paper is to present a microscopic 
theory for the phenomenological parameters which char- 
acterize (Ga,Mn)As as a ferromagnet. The form of the 
micromagnetic energy functional appropriate to the sym- 
metry of the crystal will be addressed in Sec. IV. We 
expect that these predictions will be useful in interpret- 
ing experimental studies of these semiconductor's mag- 
netic properties. The search for the carrier density, Mn 
concentration and (III,Mn)V compound semiconductor 
material for which the critical temperature has its max- 
imum value is, perhaps, the most important endeavor 
in this research p^ep-i It has been guided to date by 
mean-field-theoryE3"t3 considerations that neglect the 
low-energy correlated magnetization fluctuations charac- 
terized by Ms, r^,|-and A. However, as we have em- 
phasized earlier J2jE3 small exchange can also limit the 
temperature at wiuch. long-range magnetic order exists. 

In recent workEj~E2l we developed a theory of diluted 
magnetic semiconductor ferromagnetism which accounts 
for dynamic correlations in the ordered state. For sim- 
plicity we used parabolic bands for the itinerant carri- 
ers, leading to isotropic ferromagnetism. We found, in 
addition to the usual spin-wave mode, a continuum of 
Stoner excitations, and another collective branch ot-cp-. 
citations, optical spin waves. As we have shownj23'LHI 
the low-energy Goldstone modes can suppress the crit- 
ical temperature in comparison to mean-field estimates 
and can even change trends in Tc as a function of the 
system's parameters. 
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To address anisotropy one has to go beyond a 
parabolic-band model and take details of the band struc- 
ture into account. The itinerant-carrier bands are p-type 
and reflect the crystal symmetry of the underlying lat- 
tice. Due to spin-orbit coupling the spin degrees of free- 
dom also feel the crystal anisotropy. We will use a six- 
band description based on the Kohn-Luttinger Hamilto- 
nian to model the (Ga,Mn)As valence bands. States near 
the Fermi energy in our model include substantial mix- 
ing with the split-off valence band making the six-band 
model we epHplew-ithe minimal band model. M&an-field 
calculationfaO'Ea as well as Monte Carlo studiesEa based 
on this more realistic band structure have been performed 
recently to address magnetic anisotropy effects and ex- 
plore trends in the critical temperature. 

In Sec. ^ we set the starting point of our theory by 
deriving a formal expression for the effective action for 
the Mn impurity spins after integration out the itiner- 
ant carriesSpi This formal development generalizes ear- 
lier workcJO to general spin-orbit coupled band models. 
Using the effective action we numerically determine the 
zero-temperature spin- wave dispersion in Sec. III. Then, 
in Sec. IV we establish the connection between our spin- 



wave results and the classical micromagnetic energy func- 
tional, adjusted to the symmetry defined by the crystal 
structure. Finally, in Sec. 0, we present numerical re- 
sults for the spin-wave dispersion and the exchange con- 
stant. We find that for the parameter range of interest 
the exchange constant is enhanced by up to an order of 
magnitude compared to naive results obtained earlier, in 
which light-hole bands were neglected and the heavy-hole 
bands were approximated as parabolic. 



II. HAMILTONIAN AND EFFECTIVE ACTION 

Our theory is based on the following model. Magnetic 
ions with spin S* = 5/2 at positions R/ are antiferromag- 
netically coupled to valence-band carriers described by 
an envelope-function approach. 



H = Ho + Jpd Jd^r S(r) • s(r) 



(3) 



where S(r) — J^i S/(5(r— R/) is the impurity-spin density 
and Jpd > 0. We approximate the impurity-spin density 
by a continuous functions (instead of a sum of delta func- 
tions) . Provided that Jpd is not strong enough to localize 
valence-band carriers near the Mn sites and the average 
distance between two Mn ions is small in comparison to 
the Fermi wavelength of the itinerant carriers, this ap- 
proximation can be justified. The itinerant-carrier spin 
density is expressed in terms of carrier field operators by 
s(r) = where i and j label the basis 

in Hilbert space of spin and orbital angular momentum, 
and Sy is the corresponding representation of the spin 
operator. The envelope-function Hamiltonian Hq for the 
valence bands can be parametrized by a small number 



of symmetry- adapted parameters, the Luttinger param- 
eters 7i, 72 , 73, and the spin-orbit coupling Ago which 
splits the six states at the band edge into a quartet and a 
doublet. For explicit expressions of Hq as well as repre- 
sentations of the spin matrices in a coordinate system in 
which the x, y, and z-axis are along the crystal axes, we 
refer to Eqs. (A8)-(AlQLof Ref. |27[ based on the Kohn- 
Luttinger Hamiltonian.Lj 

Phenomenological Hamiltonians of this form have 
proven successful in understanding optical proper- 
ties of the cloaely related magnetically doped II-VI 
semiconductors. B'Q In that case however, the inclusion of 
direct antiferromagnetic interactions between Mn spins 
that lie on neighboring lattice sites proved to be im- 
portant. This interactions appear to be much weaker 
in the III-V_case, although this difference is not fully 
understood.E3 We do not include these spin-spin interac- 
tions in our calculations, because we have no knowledge 
of the strength of the coupling coefficient. If included, 
they would mainly influence short-distance spin excita- 
tions which are not in any event the focus of this work. 

We note that local-spin-depsity-approximation elec- 
tronic structure calculations j2a taken at face value, find 
a strong hybridization between p and d carriers, which is 
not completely consistent with the niodisl we use. How- 
ever, given the overwhelming successErQ of the present 
model in the case of paramagnetic (II,Mn)VI materials, 
we find it unlikely that the Mn d electrons are itinerant. 
If they were, the phenomenological approach we study 
here would be incomplete. 

Our treatment of the Mn spin system as a continuum 
greatly simplifies our calculation by eliminating disor- 
der associated with randomness in the Mn positions. 
We do anticipate that the disorder can influence both 
anisotropy and exchange constant at the low- and high- 
carrier-concentration extremes. These effects are, how- 
ever, outside the scope of the present study. We will iden- 
tify several important features of the microscopic physics 
in the disorder-free model that we expect to be robust. 

Our first goal is to integrate out the itinerant carriers 
and arrive at an effective description for the impurity- 
spin degrees of freedom. This procedure follows the 



analysis presented in Ref. 29 for the (simpler) two-band 
model. For small fluctuations around its mean-field po- 
larization, we can approximate the spin operators by 
S+{v) « 6(r)V2]VM;^, S-{r) = {S+{r))\ and S^ir) = 
NunS — bU:r)b{r) , where b^r),b{r) are bosonic Holstein- 
PrimakoflE3 fields. The quantization axis z is chosen here 
along the zero-temperature spin orientation. After inte- 
grating out the itinerant carriers the partition function, 
Z = f Vlzz] exp{—Scs[zz]), is governed by the effective 
action for the impurity spins, 

S^s[zz] = Sbp[zz] - Indet [{G^^y^ + 6G-\zz)] , (4) 

where 6*3? [zz] — Jq dr J d^rzdrZ is the usual Berry's 
phase term, and the complex numbers z and z label the 
bosonic degrees of freedom. In Eq. (Q), we have already 
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split the total kernel G ^ into a mean-field part {G^^) ^ 
and a fluctuating part 6G~^, 

(G^F)rl ^ (Or - m) 6^J |j) + A^Mn Jpd54- (5) 



5G-j\zz) = (zs-. + ZS+) ^2NunS - 2zzst 



(6) 



where /i denotes the chemical potential, and i and j range 
over a complete set of hole-band states. In the following 
we define the mean-field energy A = NunJpdS to flip the 
spin of an itinerant carrier. The physics of the itinerant 
carriers is embedded in the effective action of the mag- 
netic ions. It is responsible for the retarded and non-local 
character of the interactions between magnetic ions. 



III. INDEPENDENT SPIN- WAVE THEORY 

Independent spin-wave theory is obtained by expand- 
ing Eq. (H) up to quadratic order in z and performing 
Matsubara imaginary time and space Fourier transforms. 
Since 6G~^ is at least linear in z the series 



Scs[zz] 



(7) 



can be truncated after 7i — 2, where n denotes the or- 
der in SG^^. The zeroth-order contribution, S'^^^^ [22;] = 
SBp[zz]~'lndct{G^^)^^ contains the Berry's phase term 



Sbp [zz] 



E(- 



-Wm)z(k, t'„)z(k, 



(8) 



and the mean-field contribution from the itinerant carri- 
ers, which is independent of the bosonic fields z and z. 
Here, Vm are the bosonic Matsubara frequencies. 

The next term of the expansion, S'^^^ [zz] — 
— tr (G^^SG~^), reads in Fourier representation 



Si^^[zz] 



y^^z{k,iy„i)z(k,iy„i) , 



(9) 



m.k 



plus terms linear in z and z. Here, ujn and Vm are 
fermionic and bosonic Matsubara frequencies, respec- 
tively. To determine the mean-field Green's functions 
we diagonalize the matrix (i|i?o(q) + ^s^[j) for each 
wavevector q and denote eigenvalues and eigenstates 
by ea(q) and |a), respectively. Since the itinerant- 
carrier spin density (s) = (1/1^) Eq Eq f[<^a{q)]{a\s\a) 
is aligned antiparallel to the impurity spin-polarization 
axis (otherwise this would not be an easy axis) , the terms 
linear in z and z drop out, and we get 



213V 



E 

m.k 



z{\<i,l'm)z{\<i,Vra) ; 



(10) 



where f = —2{s^)/p is the fractional itinerant-carrier po- 
larization, < ^ < 1. 

(2) 

For the second-order term of the expansion, S^g [zz] — 
\ tr (G^FjG-iG^F(5G-i) , we find 



EE 



/[eQ(q)]-/[eMq + k)] 



4/3^^2 ^ ^ + g^(q) - e^(q k) 

m,q,k a/3 

+S-f3S'^^z{^k,lyr,^)z{~k,~l^„^) \ + O (z^) , (11) 



with s 



Q/3 



and s 



pa 



\s^\a) 



Note, that the 

indices a and (3 label the single-particle eigenstates for 
valence-band carriers with different wavevectors, namely 
q and q+k, respectively. 

The Matsubara frequency v^^ in the denominator on 
the r.h.s of Eq. ( pi] ) accounts for the dynamics of the 
itinerant carriers. This frequency dependence is crucial 
to account for the existence of the Stoner spin-flip contin- 
uum and the optical spin- wave mode. On the other hand, 
the existence of the usual spin wave follows already from 
the static limit (i.e., when the frequency dependence in 
the denominator in Eq. ([l^) is dropped), and the spin- 
wave dispersion is described rather accurately. 

The sum of Eqs. (||), (p^), and ( pT| ) is a quadratic form 
in the bosonic fields z and z. The zeros of the kernel de- 
fine the spin-wave energies ilk as a function of momentum 
k (after analytic continuation ivj^ Q + iO+). In the 
following we go the static limit as discussed above. We 
define the quantities 



El' 



q a/3 



/[eQ(q)]-/[e^(q + k)] ^ 
ea(q) - e/3(q + k) 



(12) 




with (7, ct' = ±, and perform a Bogoliubov transforma- 
tion, which eventually yields 



(13) 



From the definition Eq. (|l^) we see that is real and 
E^^ = {E^'^y . Equation ( p^ ) is the central result of 
this section. The remaining task is to evaluate the frac- 
tional itinerant-carrier polarization ^ and the quantities 
E^^ and E^^ numerically. 

Before we carry on with establishing the relation be- 
tween Eq. (|l^) and micromagnetic parameters of a clas- 
sical energy functional, we make three remarks: 

(i) Correlation effects among the Mn spins, which are 
not described by the mean-field picture, enter our theory 

(2) 

via the contribution S'^g [zz] . To reduce our theory to the 
mean- field level we would have to neglect this term, i.e., 
truncate the series Eq. already after n = 1. In this 
case, the energy — JpdP^/2 of a Mn spin excitation 
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would be dispersionless and by a factor of p^/(2iVMn'S') 
smaller than the mean-field energy A to flip an itinerant- 
carrier spin. Due to correlations between Mn and band- 
spin orientations, however, the spin-wave energy is 
always smaller than il^^ . 

(ii) In the absence of spin-orbit coupling all products 
of the form s^^s^q, or s~^s^^ vanish. As a consequence. 



= for all k. This statement is even true 



for finite spin-orbit coupling in case when the valence 
bands are isotropic. 

(iii) To go beyond the static limit, we may expand 



up to linear order 
this linear correc- 



the fraction on the r.h.s of Eq. (Ill 
in ivm- As shown in Appendix^ 
tion simply amounts to a renormalization of the Berry's 
phase term by replacing 17(1 — x). In the ab- 

sence of spin-orbit coupling we find that x is just the 
ratio of the spin densities, x — {s^)/{NMnS), where 

(...) = (l/l^)EqEa/[ea(q)](«l ■•■!«)■ For finite spin- 
orbit coupling but with a band Hamiltonian that is in- 
variant under rotation in space, the spin has to be 
replaced by the total angular momentum -I- i.e., we 
&ndx = {s'+F)/iNMnS). 

The renormalization factor {I — x) indicates that in a 
semiclassical picture, as employed in Sec. IV, the effective 
spin density is not quite given by the Mn impurities, but 
has to be reduced due to coupling to the valence-band 
carriers. On the other hand, x is always small since the 
impurity spin density Nmu is larger than the itinerant- 
carrier concentration p and the Mn spin S* = 5/2 is com- 
paratively large. We, therefore, stick to the static limit 
in the following discussion. 



IV. EASY AXIS, ENERGY GAP AND SPIN 
STIFFNESS 

The purpose of this section is to establish the connec- 
tion between the spin-wave dispersions evaluated later 
and the micromagnetic energy functional. The non-local 
magnetostatic contribution which is omitted from our 
theory can be added as needed in applications. The 
short-range part of the functional is a symmetry-adapted 
gradient expansion of the energy density e[ri]. In mag- 
netism literature the non-constant portion of the zeroth- 
order term e^"'[n] is known as the magnetic anisotropy 
energy, and the leading gradient term e'^'^fnl is known 
as the exchange energy. As shown in Ref. magnetic 
anisotropy effect in ferromagnetic semiconductors are, in 
the absence of strain, very well described by a cubic har- 
monic expansion which is truncated after sixth order, 

e'"''[n] = Ki {nlnl + n\nl + n^n^) + K2 {n^Uyn^f 

(14) 

with anisotropy parameters Ki and K2. Correlation 
of spin polarizations at different positions are described 



by the gradient term e'"^[n]. In order to address long- 
wavelength spatial fluctuations, we expand the gradient 
term up to lowest nonvanishing order. 



n = 



Aab\danb\'^ . 



(15) 



with exchange parameters Aab- We find in our numerical 
calculations that anisotropy in the exchange constant is 
negligibly small, i.e., we can choose Aab — A for all a, 6, 
as generally assumed in the magnetism literature. 

To establish the connection of the energy functional to 
our microscopic spin-wave calculation, we first have to 
determine the direction of the mean-field spin polariza- 
tion and then to study small fluctuations. The first step is 
achieved by minimizing the energy Eq. (|l4|) . It is easy to 
show that the mean-field orientation n'^''^:an only point 
along a high-symmetry axis (100), (110), (HI), or an 
equivalent direction (except for the special case K2 = 
and Ki > 0, where we find an easy-plane anisotropy in 
the planes Ux — Q, Uy — Q or — 0, and, of course, the 
isotropic case Ki — K2 = 0). In Fig. |l| wc show how the 
mean- field polarization direction depends on Ki and K2. 

Now we consider small fluctuations around the mean- 
field orientation fi^^. The kernel of the quadratic form 
e'^'"[n] — e''"'[n'^^] has two eigenvalues Ai and A2. We 
find for along (100) that Ai = A2 = Ki, for (110) 
we get Ai = -Ki and A2 = {2Ki + K2)/^, and for (111) 
we obtain Ai = A2 = —{QKi + 2i4r2)/9. In all cases, 
Ai and A2 are positive. Quantizing the collective spin 
coordinate at long wavelengths, it follows that 



2VAIA^ , 2A 



(16) 



There are two alternative ways to determine the energy 
gap 0/j=o- One can either perform the fc = limit of the 
spin-wave dispersion Eq. (^3|) or evaluate the coefficients 
Ai and A2 from calculating the energy for mean-field ori- 
entation n'^^ along the three high-symmetry axes (100), 
(110), and (111). The numerical effort of the latter pro- 
cedure, which has been used previously in Ref. |2^, is 
lesser for given accuracy. The virtue of the spin-wave 
calculation is to determine the exchange constant A. 
We conclude this section with two remarks: 

(i) In case that the easy axis is along (100) or (111) 
(as it is for all parameter sets considered in Sec. ^), the 
energy cost of tilting the polarization axis by small angles 
is independent of the direction of the deflection. This is 
required by symmetry and indicated by the fact that Ai 
equals A2 . As shown in Appendix ^ the term -Ejtjg then 
vanishes, and the fc = limit of Eq. (|l^) is identical 
to the energy gap calculated from standard perturbation 
theory where the perturbation describes the deviation of 
the spin polarization from the mean-field direction. 

(ii) The denominators NmhS in Eq. (16) corresponds to 
the static limit employed in our calculation. The renor- 
malization of the Berry's phase term due to corrections 
to the static limit (see discussion in the previous section) 
could be accounted for by multiplying NunS with (1 — a;). 
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V. NUMERICAL RESULTS FOR THE 
SPIN- WAVE DISPERSION 

In this section we present numerical results for the 
spin-wave dispersion of (Ga,Mn)As. From these calcu- 
lations we extract the spin stiffness as a function of the 
itinerant-carrier concentration p and on the exchange 
coupling Jptj. To model the sample which showed the 
highest transiticin temperature of llOK so far, we choose 
as parameterg23 -/Vmh = lnm~'^, p = 0.35 nm~'^, and 
^pd = 0.068 eVnm^'^. As a consequence, the mean- 
field spin-splitting gap for the itinerant carriers is A = 
iVM„Jpd5 = 0.17eV. 



A. Isotropic vs. six-band model 

The origin of ferromagnetism, the nature of the spin 
excitations, and trends in the critical temperatures can 
be explained within a simple model which dascribes the 
itinerant carriers by parabolic bands. ljEIIEZI For more 
quantitative statements, a more realistic description of 
the band structure should be employed. 

In Fig. ^ we show results for the isotropic model with 
two parabolic band with effective mass m* — O.Smg (a 
Debye cutoff ko with — Qii^Nyin ensures that we 
include the correct number of magnetic ion degrees of 
freedom). We find for the majority-spin Fermi energy 
tp = 0.44 eV measured from the bottom of the band. 
This yields the itinerant-carrier polarization ^ = 0.35. 
The dashed line in Fig. ^ marks the mean-field spin- 
fiip energy, obtained by neglecting correlation. Since the 
spin- wave energies are far below the mean-field result, the 
isotropic model suggest that for these parameters corre- 
lation is very important and the critical temperature is 
limited by collective fluctuations. According to the clas- 
sification scheme introduced in Ref. the system would 
be in the " RKKY collective regime" . 

By fitting a parabola at small momenta we 
find f7fc/A = Qm&Sik/kof which yields A = 
0.095 meVnm"^ = 0.015 pJm"^ 

Now we use the six-band Kohn-Luttinger Hamiltonian 
with Luttinger parameters 71 — 6.85, 72 — 2.1, and 73 = 
2.9 and spin-orbit coupling Ago — 0.34 eV. In the absence 
of the Mn ion, this model has anisotropic heavy- and 
light-hole bands with masses mt ~ 0.498me and mi « 
0.086me. The spirit of the naive parabolic band model 
is the hope that only the band with the larger density 
of states matters and the anisotropy is unimportant. We 
will see that these hopes are not fulfilled. 

We find a much lower Fermi energy ep = 0.25 eV than 
obtained for the isotropic model, and, therefore, a much 
higher itinerant-carrier spin polarization ^ = 0.73. By 
calculating the mean-field energy for a Mn spin polariza- 
tion along the three high-symmetry axes we determine 
the anisotropy parameters as Ki — 19.6 x lO^^eVnm^^ 



the easy axis is (100) and the energy gap, r2fe=o/A = 
9.2 X 10~^, is very small in comparison to the bandwidth 
of the spin- wave dispersion. 

In Fig. ^ we show the spin- wave dispersion for wavevec- 
tors k along the easy axis. We observe that the effect of 
in Eq. ( |l3| ) is negligibly small and can, therefore, be 
dropped. Furthermore, we find that the dispersion for a 
spin wave perpendicular to the easy axis cannot be distin- 
guished from the dispersion of spin waves along the easy 
axis within numerical accuracy. By fitting a parabola at 
small momenta we find fifc/A — nk=o/ A + 0.16(/c/fc£))^ 
which yields A — 2.2meVnm~^ = 0.36 pJm~^ (see in- 
set of Fig. H) . Furthermore, we see that the energy gap 
ilk=o I A determined in the previous paragraph is consis- 
tent with our spin-wave results. Employing the six-band 
model we see that for the given parameters the spin- wave 
energies are much closer to the mean-field value than for 
the iso trop ic model. According to our classification given 
in Ref. ^the system is rather in the "mean-field regime" , 
which explains the success of mean-field theory to repro- 
duce the critical temperature. 



B. Six-band model plus strain 

The lattice constants of (Ga,Mn)As and GaAs do not 
match. Since low-tcmperature molecular beam epitaxy 
(MBE) has to be used to overcome the low solubility 
of Mn in GaAs, even thick films of (Ga,Mn)As grown 
on GaAs cannot relax to their equilibrium. The lattice 
of (Ga,Mn)As is instead locked to that of the under- 
lying substrate. This induces strain which breaks the 
cubic symmetry. The infiuence of MBE growth lattice- 
matching strain on. hole bands of cubic semiconductors is 
well understood.L^Ej This effect can easily be accounted 
for adding a strain term to the Hamiltonian (we use 
Eq. (32) of Ref. 27 with strain parameters eo = —0.0028 
and r = -3.24 eV). 

We choose the growth direction to be along (001) and 
compute the energy for five different directions (100), 
(001), (110), (Oil), and (111). The lowest energy is found 
for (100), i.e., we find an easy-axis anisotropy where the 
easy axis is in the plane perpendicular tOL-the growth 
direction, in accordance with experiments. c£l From the 
expansion of the mean-field anisotropy energy for small 
fiuctuations around (100) and determining the eigenval- 
ues for small fiuctuations we estimate that the energy gap 
in the dispersion is larger than in the absence of strain by 
a factor of 1.4, i.e., still small. The spin stiffness derived 
from the curvature of the spin- wave dispersion is identi- 
cal to that in the absence of strain within the accuracy 
of our numerical calculations. We will, therefore, ignore 
the effect of strain in the following. 



and K-r — 1.6 x 10 °eVnm 



-3 



As a consequence. 
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C. Spin stiffness 

In Fig. 1^ we show the exchange constant (or spin stiff- 
ness) A as a function of the itinerant-carrier density for 
both the isotropic two-band and the full six-band model. 
We find that the spin stiffness is much larger for the six- 
band calculation than for the two-band model. Further- 
more, for the chosen range of itinerant-carrier densities 
the trend is different: in the two-band model the ex- 
change constant decreases with increasing density, while 
for the six-band description we observe an increase with 
a subsequent saturation. 

Experimental estimates for Jpd vary from 0.054 eVnm^ 
to 0.15 eV nm^, with more recent wcpk paggesting a value 
toward the lower end of this rangcEJo To address the 
dependence of the spin stiffness on Jpd, we show in Fig. ^ 
results for a two values of Jpd which differ by a factor two. 

To understand the different behavior of the two-band 
and the six-band model we recalEjO that the two-band 
model predicts A = ph^ / (8m*) for low densities p, while 
at high densities A decreases as a function of p (note 
that in Refs. ^ and ^ the spin stiffness was character- 
ized hy p — 2A). The crossover occurs near A e^?. 
The difference in the trends seen for the two- and six- 
band model in Fig. |^ is consistent with the observation 
that, at given itinerant-carrier concentration the Fermi 
energy is much smaller when the six-band model is em- 
ployed, where more bands are available for the carriers, 
than in the two-band case. Furthermore, we emphasize 
that, even in the limit of low carrier concentration, it is 
not only the (heavy-hole) mass of the lowest band which 
is important for the spin stiffness. Instead, a collective 
state in which the spins of the itinerant carriers follow 
the spatial variation of a Mn spin-wave configuration will 
involve the light-hole band, too. 



CONCLUSION 

In conclusion we present a microscopic calculation of 
micromagnetic parameters for ferromagnetic (Ga,Mn)As. 
We draw the connection of the anisotropy and exchange 
constant of a classical energy functional to the gap and 
curvature of the spin-wave dispersion. Numerical results 
for the spin stiffness as a function of itinerant-carrier con- 
centration and p — d exchange coupling are shown. We 
find that the energy gap is much smaller than the band- 
width, the spin stiffness is nearly isotropic, and strain 
does not effect the dispersion much. Furthermore, we see 
that a model with isotropic valence bands underestimates 
the spin stiffness considerably. 
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APPENDIX A: RENORMALIZATION OF 
BERRY'S PHASE TERM 

In this section, we show that the linear correction be- 
yond the static limit yields a correction to the Berry's 
phase. This accounts for the fact that the effective spin 
in a semiclassical approach does not have the length S 
of the Mn impurities by is reduced by a factor (1 — x) 
due to coupling to the itinerant-carrier spin degree of 
freedom. We expand the r.h.s of Eq. ( |Tl] ) up to linear 
order in ivm- The terms with zz and zz vanish since it is 
an odd function under the following operation: we shift 
q — > — q — k, use the fact that eQ(— q — k) = eQ,(q + k) 
and e^(— q) = e/5(q) and, eventually, exchange a ^ /?. 
Hence, we only have to deal with the term which involves 
zz. 

Assume that there is an operator Q such that 
iVM„Jpd5's+ = [H,Q] and A^Mn Jpd S's" = -[H.Q'^]- 
Then the renormalization for the Berry's phase is given 
by (1 — x) with 



1 



^^/[64q)](a|[Qt,Q]|a) (Al) 



independent of k. In the absence of spin-orbit coupling 
we choose Q = , which yields 



(A2) 



i.e., X is just the ratio of itinerant-carrier spin concen- 
tration to Mn spin density. For finite spin-orbit coupling 
the spin of the itinerant carrier is coupled to the orbital 
angular momentum. If the band Hamiltonian is invariant 
under rotation in space, we can choose Q = + Z+ and 
find 



P) 



(A3) 



In this case the correction is given by the ratio of the 
total angular momentum density of the itinerant carriers 
and the Mn concentration. 

If the valence bands are described by the Kohn- 
Luttinger Hamiltonian, however, the total angular mo- 
mentum is no longer a conserved quantity. Or, equiva- 
lently, the orbital angular momentum of the valence-band 
carriers couples to the crystal lattice, and the simple form 
Eq. (|A3|) no longer holds. 
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APPENDIX B: SPIN- WAVE DISPERSION GAP 

The goal of this section is to rederive the fc = limit of 
the spin-wave energy Eq. ( |l3| ) by standard perturbation 
theory where the perturbation describes the deviation of 
the spin polarization from the mean-field direction (100) 
or (111). 



1. Proof that = for (100) or (111) easy axis 



We start by showing that 

l^^/[ea(q)]-/[e/3(q)] 



a/3 



ea(q) - e/3(q) 



(Bl) 



in 9 does not contribute since (s^) = (s^) = 0. To ob- 
tain the quadratic order in 9 we use first-order pertur- 
bation theory for the term in Eq. ( |B3| ) and second- 
order perturbation theory for the cos ip + s'^ sin ip = 
(s+e"*"^ -I- s~e'^'-f)/2 contribution. For the former we get 

q a 

and the result for latter reads 

9^A^ |(a|s+e-*^- 



SE" = 



q a/3 



ea(q) - e/3(q) 



(B5) 



vanishes, if the mean-field polarization is along (100) or 
(111). Note that a and /3 now label the same basis states. 

Let q be a wavevector which is obtained from q 
by rotation about the z-axis (which is defined by the 
mean-field Mn spin polarization direction) by an an- 
gle (fi which respects the symmetry of the crystal. If 
the z-axis is (100) or (111), then the allowed angles are 
If e {0, 7r/2, TT, 37r/2} or {0, 27r/3, tt, 47r/3}, respectively. 
For (110) it would he (p £ {0,7r}. Due to the symme- 
try of the crystal, the spectrum at q, labeled by a (or 
/3), is identical to that at q, i.e., ~ £„. The corre- 
sponding eigenstates are connected by \a) = U\a) with 
U = exp[i(s^ -I- F)ip]. 

Since the spin operator s+ and the operator for the 
orbital angular momentum P commute, the following re- 
lation is satisfied. 



If the easy axis is along (100) or (111), then -E^t^Q = 
-EjJJg = 0, and the energy is independent of ip. The spin 
wave energy at fc = can now be obtained from the ratio 
of the energy change 6E' + 5E" and the change of the 
spin 5S = 9'^NMnS/2, which yields 



^k=a Jpd ( K 



fc=0 



(B6) 



and, as desired, we recover the fc = limit of Eq. (|l^) for 
an easy-axis direction (100) or (111). 



<^4a = {<^\U-'s+U\P){P\U-'s+U\a) 

_ -2iv + + 



(B2) 



As a consequence, the partial summation in Eq. (Bl) over 
all wavevectors q which are equivalent due to symmetry 
yields a factor (1 -|- e^*^ + er"^" + e^'^") = for the easy 
axis (100) and (1 + e-^"/^ + e-^"/^) = for (111), i.e., 
Et=o = in both cases. In contrast, the corresponding 
factor (1 -|- e~^") for (110) is nonzero, and E^^^ is, in 
general, finite. 



2. Relation between energy gap and -E^^g 

To determine the energy cost of tilting the spin polar- 
ization by small angle 9 out of the mean-field direction, 
we add to the Hamiltonian the perturbation 



H' = A 



-— -I- 9 {s^ cos(p + sint/j) 



(B3) 



and use standard perturbation theory. Here, (p G [0, 2ti) 
is the azimuth, and A = NunJpdS . The linear order 
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FIG. 1. Easy-axis direction as a function of Ki and K2. 
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FIG. 2. Spin-wave dispersion for the isotropic model 

for itinerant-carrier density p = 0.35 nm~^, impurity-spin 
concentration A^Mn = 1.0 nm~^ and exchange coupling 
Jpd = 0.068 eVnm"^ (which yields A = 0.17eV). 
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FIG. 3. Main panel: Spin-wave dispersion for the 6-band 
model for itinerant-carrier density p = 0.35 nm~"^, impu- 
rity-spin concentration Nmh = 1.0 nm~^ and exchange cou- 
pling Jpd = 0.068 eVnm"^ (which yields A = 0.17 eV). The 
momentum k is chosen to be parallel to the easy axis (100). 
Inset: Spin- wave dispersion on a log- log plot (circles) and the 
fit 9.2 X 10"^ + 0.16{k/kDf (solid line). 
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FIG. 4. Exchange constant A as a function of itiner- 
ant-carrier density p for the six-band and the two-band 
model for two diflFerent values of Jpd = 0.068 eVnm^^ (solid 
lines) and 0.136 eVnm~^ (dashed linos). The impurity-spin 
concentration is chosen as AfMn = 1.0 nm~^, which yields 
A = 0.17 eV (solid lines) and A = 0.34 eV (dashed lines), 
respectively. 
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